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Abstract. Recent developments in the theoretical evaluation of charmonium 
production in ultrarclativistic heavy-ion collisions (URHIC's) are discussed. In 
particular, the consequences of equilibrium properties of open and hidden charm 
states — accessible, e.g., in QCD lattice gauge calculations — are assessed. These 
include abundances as well as formation and dissociation rates of charmonia in 
both hadronic and quark-gluon matter. 

1. Introduction 

Hadrons containing charm quarks play a special role in ultrarelativistic heavy-ion 
physics [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]. On the one hand, the charm (c) quark carries a 
mass, m c ~1.3 GeV/c 2 , large enough for its production to be dominated by primordial 
nucleon-nucleon (N-N) collisions [12]; on the other hand, m c is small enough to allow 
for significant interactions with surrounding light quarks and gluons, thus providing 
a rather direct probe of the hot and/or dense matter created in central collisions 
of heavy nuclei (A- A). This particularly applies for c-c bound states (charmonia): 
based on the production systematics in p-p reactions, where they amount to ^2% of 
the total charm (cc) yield, charmonia have been predicted to be suppressed beyond 
standard nuclear effects if a deconfined medium is formed in A- A collisions [2] . At the 
same time, reinteractions of c-quarks imply approach towards thermalization, and, in 
principle, the presence of the backward channel in charmonium dissociation reactions, 

Xx + J/ip^X 2 + c + c {D + D) . (1) 

To identify the relevant processes in a heavy- ion reaction, several questions arise, e.g., 

• What are the inelastic charmonium cross sections, and corresponding equilibra- 
tion times, in the Quark-Gluon Plasma (QGP) and hadronic phase? 

• What are the in-medium properties (mass and width) of the charmonium and 
open charm states in both phases? 

• How does the primordially produced number of charmonium states compare to 
thermal equilibrium abundances? 

• What is the significance of the critical temperature T c in the previous items? 

§ email: rapp@nordita.dk 

|| email: loicOtonic. physics.sunysb.edu 



Theoretical Overview on (Hidden) Charm in High-Energy Heavy-Ion Collisions 2 



• Are there large corrections to the expected open charm content of the system, 
and do c-quarks (D-mesons) thermalize? 

In this talk we will address (some of) these issues as follows. In Sec. 2 we briefly 
recall the basic features of open charm production, including new insights from recent 
RHIC data. In Sec. 3, we first assess in-medium properties of open and hidden 
charm states from recent lattice QCD results. Second, we elaborate on their possible 
consequences on the description of charmonium evolution in URIHC's in both QGP 
and hadronic phases. The main tool will be rate equations with their main ingredients 
- initial conditions, relaxation times and equilibrium limits - being discussed in some 
detail. Third, we will comment on statistical hadronization models, in particular how 
in-medium modifications of charmed hadrons may affect their results. In Sec. 4, we 
confront various approaches to available data, and in Sec. 5 we conclude. 

2. Open Charm Production 

The production of cc pairs in heavy-ion reactions is expected to behave as a hard 
process, i.e., scale with the number of primordial N-N collisions. While p-A data 
follow this expectation [13, 14], only indirect measurements via semilcptonic decays, 
c(c) — > l + (l~)X (within charmed hadrons), are available for the A- A case so far. 

In Pb-Pb at SPS energy, ^=17.3 AGeV, NA50 [14] found an enhancement of 
intermediate-mass (IM, M w =1.5-2.5 GeV) dimuon pairs over the expected Drell-Yan 
and open charm sources, gradually increasing with centrality up to a factor of ~2.5 
in most central events. This excess is well reproduced if an increase in open charm 
production by a factor of ^3.5 is postulated. However, underlying mechanisms for 
such an increase are not easily conceivable. In fact, it has been found [15, 16, 17] 
that thermal dimuon radiation from a fireball with reasonable initial temperatures, 
Ti=200-250 MeV, can, maybe more naturally, account for the IM dimuon excess. The 
experimental resolution of this problem is expected from the NA60 experiment [18]. 

In Au-Au at RHIC energy, y / s=130,200 AGeV, single-electron (e ± ) transverse- 
momentum (p t ) spectra have been measured at various centralities [19, 20]. After 
subtraction of light-hadron decay sources, the remaining spectra can be accounted 
for by "standard" charm-production extrapolated from N-N event generators. On 
the one hand, this limits the possibility for an appreciable charm enhancement. On 
the other hand, it is somewhat surprising that even for central collisions the spectral 
shape is reasonably well reproduced without any reinteractions, as one would naively 
expect a softening if charm quarks (partially) thermalize. However, thermalization 
also implies that c-quarks participate in the collective matter expansion. In Ref. [21] 
it has been shown that the current PHENIX data [19, 20] for "charm-like" e ± arc 
also consistent with the assumption of complete c-quark thermalization and collective 
flow. As discussed below, c-quark reinteractions have important consequences for 
charmonium production; the obvious observable to disentangle the two extremes is 
the c-quark elliptic flow [21]. 

3. Charmonium in QCD Matter 

3.1. Lattice QCD Results 



Lattice QCD calculations for the finite- T free energy of a heavy-quark pair [22], 

Fqq{v,T) = Vqq (r, T) — TSqq (r, T) , (2) 
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Figure 1. Left panel: free energy of a heavy-quark pair from unquenched 
(iV/=2+l) lattice calculations [22] as a function of Q-Q separation r for various 
temperatures T. Right panel: quenched lattice results [24] for the r] c spectral 
function above and below T c . 



including dynamical fermions, are displayed in the left panel of Fig. 1. The decrease 
of the asymptotic plateau value, Fqq(T, r — > oo), can be attributed to an in-medium 
reduction of the open charm threshold in the heavy-quark potential, Vqq, and has 
been used to infer "dissociation" temperatures Tdi SS <T c for ip' and \ states, as well 
as Tdiss— 1-1 T c for the J/ip [23]. It should be emphasized that the temperature 
dependence of the plateau is remarkably continuous, even across T c ; this has interesting 
implications which we will return to below. If bound states move above the threshold 
for spontaneous decay, it does, however, not necessarily mean that they disappear 
from the spectrum, but they can rather survive as resonances^]. This is indeed found 
in two recent (quenched) lattice calculations [24, 25], where reconstructed r] c and 
J/tp spectral functions exhibit well-defined states up to T~1.5 T c , see right panel of 
Fig. 1. The masses appear to be essentially T- independent (quite contrary to the open 
charm states), whereas the widths, T nc ^, increase significantly, to about 0.2-0.5 GeV 
slightly above T c [24, 25]. This is not much different from lifetimes obtained using 
gluon-induced J/ip dissociation cross sections [26] at comparable temperatures. 



3.2. Rate Equations in the Quark- Gluon Plasma 

Let us now turn to the description of the time evolution of J / tp abundances relevant 
for heavy-ion reactions, first focusing on the QGP phase. A corresponding rate 
equation for the number of J/ip mesons in a system (fireball) of volume Vfb, 
N^{t) — ti^(T)Vfb(t) (n^(T): J/ip density), is of the form 

± Nf = -N f L(r) + G(t) = -- K(t) - N?(rj\ , (3) 

which applies as long as a well-defined J/ip state exists, i.e., t^ 1 (T)—T^(T)<^m^ . 
The second equality in Eq. (3) follows under the simplifying assumption that the 
surrounding partons (including the charm quarks) are in thermal equilibrium. The 

If note that, in a heat bath, the distinction between bound state and resonance becomes immaterial; 
all information (mass and width) is encoded in the spectral function. 
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Figure 2. Left panel: Equilibrium J/ip abundances (assuming my, (T) = m^ ac ) in 
an isotropic, adiabatically expanding Quark-Gluon Plasma at fixed (anti-) charm- 
quark number, N C c, resembling conditions at RHIC and SPS for two values of the 
charm-quark mass. The arrows mark the numbers, N^ rlm , of J/ip's produced in 
primordial (hard) N-N collisions (as extrapolated from elementary p-p reactions). 
Right panel: J/ip equilibrium numbers in the hadronic phase (under equivalent 
conditions as in the left panel) using free and in- medium charmed hadron masses. 



charm-quark and J /ip equilibrium densities are then given by 

n c (T, 7c) = n- c (T, 7c ) = 6 7c J ^f C ( m c, T) (4) 

r#(T, 7 c) = 3 7c 2 | J^Lf^ni^T) . (5) 

Here we have already incorporated the possibility of chemical q/J-equilibrium by 
introducing the charm-quark fugacity 7c . In practice, the latter will be matched to 
the total number, iV cg , of cc pairs produced in hard (primordial) N-N collisions [27], 

N c5 = \N op ^4 + Vfb ^ >( T ) 2 E ( 6 ) 

with N op =2 Vfb(T) 7c(T) n c (T). Fig. 2 (left panel) shows equilibrium J ftp numbers in 
the QGP under conditions expected for central Pb-Pb/ Au-Au collisions at SPS/RHIC. 
One observes a large sensitivity to the (in-medium) charm-quark mass, with larger 
values for m c requiring a larger fugacity 7 C (at fixed N cS ), and thus entailing a larger 
equilibrium abundance of charmonium states. 

Besides the equilibrium J/ip number the other essential quantity figuring into 
Eq. (3) is the (chemical) relaxation time (related to processes given in Eq. (1)). In 
a perturbative framework, it is determined in terms of parton-induced J /if) breakup 
(inelastic) cross sections, <r™ e ', via convolution over thermal parton distributions /*, 

The use of the standard gluodissociation cross section [1, 28], depicted in the upper left 
panel of Fig. 3, leads to relaxation times t^>10 fm/c for temperatures T<200 MeV, 
steeply falling to below 1 fm/c at T>300 MeV. If, however, according to the lattice 
data in Fig. 1, the threshold energy of the outgoing open charm states is substantially 
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Figure 3. Left panel: parton-induced J /if) dissociation processes (upper panel: 
the equivalent of photodissociation induced by gluons; lower panel: quasifree 
destruction); Right panel: Corresponding J/rjj relaxation times. 

reduced and approaches the J/ip mass, this process becomes kinematically unfavo- 
rable, see dotted line in right panel of Fig. 3. Under these circumstances, "quasifree" 
dissociation [26] , which is characterized by an extra outgoing parton (lower left panel 
in Fig. 3), is much more efficient, cf. solid line in the right panel of Fig. 3. 

Nonperturbative approaches to the dissociation of J/^'s in a QGP often involve 
the notion that also "D"-meson states survive the phase transition [8, 29]; e.g., 
in constituent-quark based pictures [29], D-mesons exhibit the "Mott" effect, i.e., 
become broad resonances in the c-q continuum above T c , thus opening phase space 
for J/ip -> "L>" + "L>" decays. 

The final ingredient needed for an actual solution of the rate equation is the 
initial condition, N®, usually taken from hard production systematics (modified by 
initial-state nuclear absorption to account for " pre-equilibrium" dynamics). One can 
then specify under which conditions the gain term in Eq. (3) can be neglected in 
applications to heavy-ion reactions, namely: either (i) 3> tqqp {tqqp'- lifetime 
of the QGP phase) and iV° ^> N^ q or (ii) no J/tfj state is supported in the QGP. 
The arrow markers in the left panel of Fig. 2 indicate how the expected primordial 
production numbers of J/^'s compare to the thermal equilibrium levels. Whereas at 
SPS energies N^ q <C N^ rim , this no longer holds for RHIC energies, i.e., a neglection 
of the gain term appears to be justified at SPS, but this is not obvious at RHIC. 

3.3. Hadronic Matter: Statistical Production and in-Medium Effects 

Based on the success of thermal models to describe the production of hadrons 
containing u, d and s quarks [30], it has been suggested [27] that also charmed hadrons 
form by statistical hadronization (of pre-existing c, c quarks) according to thermal 
weights at T c . This, in particular, has been applied to charmonium states [27, 31, 32] 
by employing Eqs. (5) and (6) but with the fugacity 7 C determined by the hadronic 
open charm states, N op — Vfb(T c ) j c ^2 a n a (T, /j,b), where a= D, D, D* , D* , A c , 
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Figure 4. n J/ip — > DD,D* D,DD* dissociation cross sections in various 
approaches; band: QCD sum rules [34], dotted line: short-distance QCD [42], 
dashcd-dotted line: meson-exchange [40], dashed line: nonrelativistic constituent 
quark model [36]. 

. . . , runs over all known open charm hadrons, and Vfb(T c ) is the hadronic fireball 
volume after completion of the phase transition (see also Ref. [33] ) . 

In the language of the rate Eq. (3), the statistical hadronization approach cor- 
responds to having reached complete thermal and (relative) chemical equilibrium at 
T c . Interactions in the subsequent hadronic phase may still alter the charmonium 
abundances, as determined by inelastic hadronic reactions of type 7r J /%p <-> DD* +cc, 
pJ/^p <-> DD, etc., which have been intensely investigated, mostly with focus on 
the dissociation direction. An example of the current status is shown in Fig. 4, 
taken from Ref. [34]: the band, which is the main result of that work, represents 
the pion-induced dissociation cross-section as assessed within a QCD sum rule 
analysis. It is in reasonable agreement with evaluations from both constituent- 
quark [35, 36] and meson-exchange models [37, 38, 39, 40, 41] (once soft form factors 
with cutoffs A<1 GeV are included), but significantly above earlier short-distance 
QCD calculations [42] based on gluodissociation [28] folded with the (very soft) gluon 
distribution functions within light hadrons. p- [36, 38] and nucleon-induced [43] 
break-up cross sections tend to be somewhat larger. But even under rather extreme 
assumptions, e.g., an average cross section of 1 mb and a total hadron density 
?i/iad=3no— 0.5 fm -3 , the chemical relaxation time, T^,=[(crt')n/ ta£ ;] _1 ~30-40 fm/c, is 
well below typical hadronic fireball lifetimes. 

The situation may change if additional medium effects are present, most notably 
reduced D-meson masses, as indicated by the lattice results in the left panel of Fig. 1, 
QCD sum rule estimates [44] or model calculations [45]. In a constituent quark picture, 
such modifications are attributed to the light quarks losing their mass, thus directly 
reflecting on (partial) chiral symmetry restoration. From a practical point of view, the 
lowering of the DD threshold has several consequences, e.g.: (i) excited charmonia 
may move above threshold (well) below T c [23] , (ii) increased phase space accelerates 
the inelastic reactions (although wave function nodes in excited charmonia can lead 
to strong deviations from naive extrapolations [46]) (iii) the equilibrium abundance of 
charmonia in hadronic matter decreases (at fixed N cS ) . The latter effect is illustrated 
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in the right panel of Fig. 2, and is completely analogous to varying m c in the QGP, 
see discussion after Eq. (6). Even a moderate reduction of the D-meson masses by 
80 MeV at T c =170 MeV lowers the J/ip equilibrium level by a factor of 2. The other 
remarkable feature of the hadronic equilibrium J/ip numbers is that they increase 
with decreasing temperature - again a consequence of the increasing 7 C at fixed N cS 
together with < 2niD + - For heavy- ion reactions this implies that, should the J/ip 
number reach equilibrium at or around T c , subsequent hadronic reactions would favor 
additional formation over dissociation (this could be even more relevant for excited 
charmonium states which only exist in hadronic matter). Such a behavior has indeed 
been observed in kinetic [47] and transport [48, 49] model calculations. 

34. Continuity? 

Before we turn to quantitative model comparisons to heavy-ion data, let us reiterate 
ramifications of the lattice QCD results [22] for the heavy-quark free energy (left 
panel of Fig. 1) on the J/ip evolution within the rate-equation framework. As 
mentioned above, the asymptotic plateau, F cS (T, r — > 00), suggests a gradual decrease 
of the open charm threshold with increasing temperature [23] , which is surprisingly 
continuous even when passing through the phase transition. At the same time, the 
J/ip mass is essentially constant. The latter implies that, in thermal and chemical 
equilibrium (j c =l), the J/ip density, Eq. (5), only depends on temperature and 
thus carries no notion of phase transition dynamics. Of course, in a finite system, 
the latent heat released in going from quark-gluon to hadronic degrees of freedom 
would imply the total J/ip number to increase according to the volume expansion 
at T c . However, this is not the situation envisaged in a heavy-ion reaction, where, 
in addition, one expects the total charm-quark number, A cE , to be approximately 
constant, encoded in the chemical off-equilibrium parameter 7 C ^1. Hence, through 
Eqs. (5) and (6), J/ip abundances are coupled to N cS and spectral distributions of 
open-charm states (figuring into the densities), parametrically as N$ oc ^/(Vps^op) 
with n op ~2n c {m* c ,T) or J2 a =D d* n a{ m * a ^T), corresponding to QGP or hadron gas 
(HG), respectively. It furthermore seems conceivable that the open-charm spectrum 
- not only the lowest level - changes rather continuously through T c , and therefore 
n^ GP ~n^, G [e.g., with in- medium charm-quark masses m*>m c =1.15-1.35GeV/c 2 due 
to heavy-quark thermal correlation energies in the QGP). In this case, the equilibrium 
charmonium numbers decrease during the phase transition according to the volume 
increase at T c , quite contrary to the chemical-equilibrium scenario (j c =l). 

4. Model Approaches at SPS and RHIC 

At SPS energies, the interpretation of the NA50 data [50, 51] is still under debate. 
On the one hand, "comover" models [52] with a constant dissociation cross section of 
cr*" ei ~l mb are able to explain the observed suppression. One should note, however, 
that the underlying comover densities reached in the early stages of a collision - 
albeit not necessarily thermal - correspond to energy densities well above critical 
ones obtained from lattice QCD. On the other hand, the recent idea [27, 31] of 
statistical coalescence of c and c quarks at T c has stimulated further activity. One of 
its strongholds is the ip' /ip ratio [27], which experimentally [53] tends to saturate at 
the thermal value in sufficiently central A- A collisions. To describe the J/ip data by 

+ In the academic limit T — > 0, all cc pairs would be in the lowest possible state, i.e., rj c mesons! 
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Figure 5. Recent model approaches to NA50 J/ip data [50] from Pb-Pb collisions 
at SPS (-y/sjvjv=17.3 GeV). Left panel: statistical coalescence [32]; right panel: 
two-component model [26]. 



statistical production alone (assuming all primordial J/i/>'s to be suppressed or not to 
form), requires an open charm enhancement by a factor of ~3 (in line with the excess 
inferred from IM dimuon spectra, cf. the discussion in Sec. 2), see left panel of Fig. 5. 
In Ref. [26] a two-component approach, combining suppression mechanisms (nuclear, 
QGP and hadronic) with statistical production at T c , has been suggested which allows 
to incorporate peripheral A- A and p-A collisions and refrains from introducing an 
anomalous charm enhancement. In this case, QGP suppression is identified as the key 
feature at SPS, see right panel of Fig. 5. 

Let us finally confront predictions for RHIC energies with (preliminary) data 
from PHENIX [54]. The left panel of Fig. 6 shows the statistical model yields [27, 32], 
illustrating inherent uncertainties due to the thus far not accurately known open- 
charm cross section. In the right panel of Fig. 6, the results of the two-component 
model [26, 55] indicate that for (semi-) central collisions at RHIC the thermal 
production at T c indeed dominates over the primordial contribution, as the latter 
is subject to strong QGP suppression [11]. This has interesting consequences for the 
excitation function from SPS to RHIC energies [26]. The upper solid line in the right 
panel of Fig. 6 is obtained from calculations [56] solving the full rate Eq. (3), using the 
gluodissociation cross section with a vacuum J/-0 mass and width. Although is 
probably increased, the idea of a J/ip state in the QGP is supported by lattice QCD, 
recall the right panel of Fig. 1. The J/ip formation in this calculation (which does not 
invoke a hadronic phase) indeed occurs preferentially in the early QGP stages which 
accounts for most of the difference to the 2-component model (where, in turn, J/ip 
formation in the QGP has been neglected). Predictions in line with the current data 
were also made within a multiphase transport model in Ref. [48], see also Ref. [49] 
for a recent analysis. In both these calculations, secondary J/tp formation, i.e., the 
backward reaction in Eq. (1), again plays a decisive role. In Ref. [57] it has been 
argued, that a reduced nuclear absorption could be consistent with the data without 
invoking any regeneration processes. 
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Figure 6. Model predictions for J/tp production at RHIC (y / S]v r iv=200 GeV) 
compared to preliminary PHENIX data [54]. Right panel: statistical hadroni- 
zation [32]; left panel: dynamical approaches (lower solid line: two-component 
model [26], upper solid line: yields from solving rate equations in a QGP [56]). 



5. Conclusion 

About 20 years after the first ideas on utilizing charmonium states as a probe of 
hot and dense QCD matter in heavy-ion collisions, important insights continue to 
emerge. In addition to experimental data, lattice QCD calculations are now providing 
valuable information for theoretical model approaches, most notably the survival of 
(some) charmonia as resonance states in the QGP as encoded in pertinent spectral 
functions, and the smooth behavior of open charm thresholds through both hadronic 
and early QGP phases, across T c . We attempted to infer consequences thereof in the 
heavy-ion environment, e.g., the necessity to include charmonium formation reactions, 
the behavior of J /ip equilibrium levels across the hadronization transition and its 
interplay with in-medium effects on open charm states. The current high-energy 
frontier at RHIC (and future LHC) implies, for the first time, multiple production 
of cc pairs, as well as more favorable conditions for their reinteractions. If charm 
quarks indeed thermalize in the course of an A- A collision, the equilibrium level of 
charmonium states being comparable to primordial production, necessarily entails 
significant regeneration. Whether such a regime has already been reached with 
RHIC energies cannot be decided from the present J /if) and single-e 1 * 1 data. Future 
measurements in the upcoming RHIC runs (including open charm observables and 
possibly excited charmonia), in connection with further theoretical developments, will 
hopefully hold (some) answers. 
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